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We propose a graphene-based model for realizing a new type of gapless condensate by pairing 
of electron-like (n) carriers of a Dirac cone conduction band with hole-like (p) carriers of a Dirac 
valance band. Ferromagnetic superconductivity (FS) in monolayer graphene or pairing between 
oppositely (n and p) doped layers of a double layer graphene allow for the formation of this p-n 
superconductivity. For FS in graphene, the p-n condensate dominates the zero temperature phase 
diagram at low levels of doping and high exchange fields. We show that p-n pairing with p+ip-wave 
symmetry presents a stable condensate phase, which can cover the phase diagram up to surprisingly 
strong exchange fields. Our study reveals that the characteristics of relativistic quantum physics 
affect the interplay between ferromagnetic ordering and superconductivity in a fundamental way. 

PACS numbers: 74.78.-w, 12.38.-t, 75.75.-c, 74.70.Wz 



The exclusive phenomena of superconductivity and fer- 
romagnctism arc two of the fundamental macroscopic 
manifestations of quantum physics governing the state 
of electrons in a conducting material. While the former 
is based on the binding of electrons with opposite spin 
into the so-called Cooper pairs, the later tends to align 
spin of electrons in order to suppress Coulomb repulsion 
through Pauli exclusion. The old and long-standing^ in- 
terest in their interplay has been greatly intensified dur- 
ing the past two decades. The motivation has come from 
the experimental realization of artificial nano-scalcd hy- 
brid structures of ferromagnets and superconductors — , 
which allow for a controlled study of interplay between 
these competing correlations, and also by the discovery of 
new magnetic superconducting materials, iron-based su- 
perconductors^ being the most recent ones that exhibit 
high transition temperatures and exciting properties. 

Until now, most of the studies on ferromagnetic su- 
perconductivity (FS) focus on Bardccn-Coopcr-Schricffcr 
(BCS) limit — ~— in ordinary electronic material with a 
high chemical potential, in which the electrons, as non- 
relativistic particles, occupy essentially parabolic energy 
bands having an almost constant density of states in the 
scale of the superconducting pairing gap. In this respect, 
the question of FS for relativistic fcrmions has remained 
unanswered. In addition to its emergence in the context 
of high temperature and heavy fermion superconductiv- 
ity — , color superconductivity in dense quark matter at 
low temperatures^ and neutron-proton pairing in nuclear 
matter—, the question has become of particular impor- 
tance and relevance by the recent discovery of graphene 
11 i 12 , in which electrons at low energies behave like mass- 
less Dirac fermions with a pseudo-relativistic chiral prop- 
erty. One anticipates that the features of relativistic 
quantum physics to affect the magnetic ordering^ and 
the superconducting pairing^ in a fundamental way, and 
hence their mutual interplay and evolution. The aim of 
the present study is to address this question within a 
model based on graphene. We present zero temperature 
phase diagrams of relativistic FS in graphene with s-wave 
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FIG. 1: (color online) (a) Schematic of p-n pairing (illustrated 
by wavy lines) between electron-like carriers of spin a in the 
Dirac cone conduction band of the if- valley and hole- like car- 
riers of opposite spin in the valance band of the K' -valley. The 
pairing can be realized in a ferromagnetic monolayer graphene 
(b, upper panel) with an exchange field h exceeding its chem- 
ical potential fx, or between layers of a double layer graphene 
(b, lower panel) with a difference in the chemical potentials 
Afi (can be controlled by a perpendicular electric field E ) of 
the two layers exceeding their mean value (fii + fi2)/2. The 
zero-temperature phase diagram of graphene-based ferromag- 
netic superconductivity showing dependence on h and fi for 
s-wave (c) and p+ip-wave (d) pairing symmetry. The super- 
conducting region consists of three phases of BCS, Sarma and 
p-n pairing. The dashed line in (c) and (d) indicates the insta- 
bility boundary for BCS phase separating stable BCS states 
(left side) with the unstable (right side) ones. The S to N 
transition on these lines is first order. While the solid line in 
(d) indicates a second order S-N transition. 



and p+ip-wave symmetry (sec Fig. [T};,d) which explore 
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their fundamental difference from those of an ordinary 
FS. We show that FS in graphene provides the possibil- 
ity for realizing a new gapless condensate phase in which 
Cooper pairing occurs between chiral carriers of oppo- 
site types, electron- like (n) and hole-like (p), which be- 
long to different Dirac conduction and valance bands (see 

Fig. at). 

In the case of an ordinary s-wave superconductor, the 
presence of an spin-splitting field — £ causes a mismatch- 
ing of the Fermi surfaces of the spin-down and spin-up 
electrons, and hence suppresses the number of paired 
states. This can lead to a first order quantum phase 
transition from superconducting (S) state to the normal 
(N) state in the so-called Clogston-Chandrasckhar limit 
specified by a critical field, h c = Ao/y/2, with Ao be- 
ing the S pairing gap at zero temperature and zero field. 
For the fields below h Ci FS exhibits properties of a ho- 
mogeneous BCS condensate. There can also be existed 
homogenous S state with A s < h, known as Sarma or 
breached pairing state — , which involves pairing between 
the spin subband with a smaller Fermi surface and that 
with a larger Fermi sphere by curving out of a gap in- 
side the larger Fermi sphere. This leaves the exterior 
part of the larger Fermi sphere unpaired, and give rise 
to the gapless property. Sarma state exhibits an insta- 
bility, which can be avoided, for instance, by allowing a 
finite-range momentum dependence of attractive inter- 
action between fermions in a single band Fermi system 
— , and interband pairing in the case of two bands Fermi 
system — . Here we demonstrate that in graphenc-based 
FS, the new relativistic p-n superconducting phase dom- 
inates over Sarma state for s-wave symmetry, and over 
both BCS and Sarma states for p+ip-wave symmetry at 
low levels of doping. We find that, unlike the s-wave 
p-n phase, the corresponding p+ip condensate is stable 
and can be preserved in much strong exchange fields as 
compared to a conventional FS state. 

To be more specific, we consider a graphene sheet in 
the presence of a spin-splitting exchange field h and the 
superconducting pairing interactions with an on-site cou- 
pling constant go and off-site nearest neighbor coupling 
constant g±. In addition to the possibility of an in- 
trinsic superconductivity, with plasmon or phonon me- 
diated pairing interactions, in graphene—, progress have 
already been made in proximity-inducing superconduc- 
tivity by fabrication of transparent contacts between a 
graphene monolayer and a superconductor (see for in- 
stance, Refs^). Similarly, ferromagnetic correlations 
can be induced by using an insulating ferromagnetic sub- 
strate or by F metals or added magnetic impurities on top 
of the graphene sheet — . The induced magnetization is 
expected to be directed in-plane of graphene sheet, which 
allows for neglecting the orbital effects^. 



we consider the following tight-binding Hamiltonia: 
H = - y^(/i + (rh)n la - t (atybjer + h.c.) 
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where the operators ai a ja i iG ibi a /b ia ) are on-site cre- 
ation/annihilation operators for electrons of sublatticc 
A (B) with spin a — ±, respectively, t ~ 2.8 eV is the 
hopping energy between nearest neighbor carbon atoms, 
hia- is the on-site particle density operator, and fj, is the 
graphene chemical potential (we use the units such that 
fi = 1 = fcs). We constrain ourselves to the case of zero 
temperature, where the mean-field approximation in two 
dimensions is expected to have the highest validity^. In- 
terestingly, we have noticed that the Hamiltonian (JTJ) can 
equivalently describe BCS pairing of electrons from dif- 
ferent layers of a double layer graphene, with levels of 
doping jUi and \i2 in the layers 1 and 2 (see Fig.[T]D, lower 
panel). In this respect, the difference A/i = (^ — 1^2) 
and the mean (fj,\ +fj,2)/2 values in the double layer prob- 
lem play the same rule as the exchange field h and the 
chemical potential /1 in FS, respectively. Thus, interlayer 
superconductivity in double layer graphene will exhibit 
the same physics as FS in monolayer graphene. 

The on-site and the nearest neighbor parts of the in- 
teraction in (fTJ can be decoupled by the following spin 
singlet superconducting order parameters which have, re- 
spectively, s-wave and p + zp-wave symmetries 
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Replacing these mean field order parameters in H and 
performing a Fourier transformation, we obtain the 
Hamiltonian in the space of 2D wave vector k of elec- 
trons. The excitation spectrum can then be obtained by 
a Bogoliubov diagonalization as 



Ki = ^J(^ + ^) 2 + (A s + ^A p /ty-ah, (4) 

where I = ±1 denote two branches of the spectrum and 
£k+K<'> = ^pk, in the linear approximation of the free 
electron spectrum around the two Dirac points K and 
K'; vp = y is the Fermi velocity of Dirac particles, and 
k is the magnitude of k. 

From the excitation spectrum, we calculate the ther- 
modynamic potential f2g using the relation 



Slf. 



E 



rf 2 k 

(2^ k; 



2e(-££) - 1 



(5) 



To describe the on-site and off-site FS of 7r electrons in 
the honeycomb lattice structure of a monolayer graphene, 



where Eq = A^/g + 3Ap/gi is the condensation energy 
and Q(x) is Hcaviside Theta function. In the following 
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FIG. 2: (color online) The plot of quasi-particle excitation 
spectrum _E^ ; versus in the superconducting state. The 
dispersions correspond to the regimes of (a) BCS pairing with 
no effective Fermi surface for both branches I = ±1, (b) Sarma 
pairing with two effective Fermi surfaces for the branch / = — I 
but no Fermi surface for the other branch 1=1, and (c) p-n 
pairing with single Fermi surface in each branch / = ±1. 



we will see that the presence of O term in Eq.©, which 
depends on whether E^ is positive or negative, has cru- 
cial effect on the type of the pairing resulting from the 
gap equation. The conditions that cause E£ t to be ei- 
ther positive or negative can be expressed in terms of the 
topology of the effective Fermi surfaces, as we investigate 
now. From Eqs.(U]), one can see that for h > 0, A s > 
and A p > 0, and for any value of I, and \i the spectrum 
E kl is always positive. This situation is similar to that of 
a conventional superconductor in which exciting a paired 
state requires a finite energy. In contrast, the spectrum 
E kl is not positive always. Depending on the involved 
parameters, it can undergo change of sign once for I = 1 
and either once or twice for I = — 1. In Fig. [2]the possible 
excitation spectrums are shown. By analyzing the roots 
of the equation Elr = 0, we find three distinct regimes 
for h with their own Fermi surface topologies. 
In the range 



h < 



fxA p /t\ 



v/l + (A P A) 2 

there is no effective Fermi surface for both I = ±1, as can 
be seen in Fig. [2ji. In this regime the pairing is of BCS 
type for both s-wave and p+ip-wave cases. For A p = 
and A s ^ 0, this regime corresponds to an isotropic BCS 
superconductivity. Including a nonzero A p , can decrease 
slightly the area of this BCS phase in the phase diagram 
of s-wave FS. 

The second regime occurs in the range 

l^ML < h < v^ta!, 

where the pairing is of breach-pairing type. For 1=1, 
there is no effective Fermi surface as in BCS regime. 
While for the branch I = — 1 there are two roots for 
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FIG. 3: (color online) (a) Nonzero solutions of s-wave pair- 
ing gap equation versus exchange field h for different levels of 
doping /j,. We have set A s o = A/100, (b) Difference in ther- 
modynamic potentials of N and S states versus h calculated 
at the same values of as in (a). 



El, = 0, and correspondingly two effective Fermi sur- 
faces as is shown in Fig. ^jp. This regime represents to 
Sarma phase in the phase diagram. For A s = and 
A p 7^ 0, as yti — > 0, the corresponding region in the phase 
diagram vanishes. Also for A s ^ and A p = 0, this 
Sarma phase vanishes as [i — > (sec Fig. QJi,b). 

The third range h > a/^ 2 + A^ characterizes a new 
regime for which each of / = ±1 branches contains one 
effective Fermi surface (see Fig. [2}:). Importantly, one 
can see that the superconductivity in this regime is re- 
sulted by pairing between normal Fermi surfaces of carri- 
ers of opposite spin, when they lie in different type, con- 
duction or valence, subbands. In this regime ferromag- 
netic graphene presents a spin chiral material, in which 
opposite-spin carriers are of different types, electron-like 
and hole-like^. 

According to the above three ranges for exchange field, 
the gap equation has different solutions, which we analyze 
them in the following. We obtain the gap equations by 
minimizing the thermodynamic potential Eq. ([S]) with re- 
spect to the order parameters A s and A p , dfl/dA SiP = 0. 
In general, these are two coupled equations, from which 
the self-consistent solutions of the gaps can be calculated. 
We investigate the s-wave and p+ip-wave solutions sepa- 
rately, so in this case gap equations are not coupled. Let 
us start with the case when only the on-site attractive 
interaction exists (go ^ 0) and g± — 0. In this case, for 
s-wave condensate from Eqs.([5]) we obtain the gap equa- 
tion as 
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Similarly for the off-site nearest neighbor superconduc- 
tivity (go = 0, g\ 7^ 0), with p + ip-wave symmetry, we 
obtain 




d 2 k g QjEl) N 
(2tt)2 3 El , + h. 



A p = 0. (7) 



In Fig. [2K the nonzero solutions of Eq. ^ as a func- 
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FIG. 4: (color online) The same as in Fig. [3l but for p+ip- 
wave symmetry. Here we have set A p o = A/100. 



tion of h are shown for different levels of doping. At high 
dopings /i 3> Aso, the situation is similar to an ordinary 
s-wave FS wherein only BCS and Sarma solutions are 
possible. The BCS solutions extent over < h < A s o 
with a constant gap A s = A s o- The Sarma solutions 
appear in the range 1/2 < h/A s o < 1, for which the 
pairing gap increases with increasing h. Lowering \x has 
no effect on BCS solutions, but leads to a su ppression of 
Sarma states. At a critical doping with fi = \J h 2 — A 2 p- 
n pairing solutions also begin to appear within a sharply 
narrow range of the exchange field. Lowering fx further 
leads to a growing of p-n states and further suppress- 
ing of Sarma solutions. Approaching Dirac point /i = 0, 
p-n states dominates fully over Sarma states. We have 
analyzed the stability of these solutions in Fig. [3b, by 
plotting the difference in thermodynamic potentials of N 
and S states, Qs~n, as a function of h for different val- 
ues of fj,. We see that only BCS solutions arc stable and 
that is up to the critical field h c , wherein graphene FS 
undergoes a first order phase transition into N state. In- 
terestingly, the range of stability of BCS states increases 
with decreasing the doping and tends to be over the full 
range of < h < A s o in the limit of /j, <§; A s o- We note 
that Sarma and p-n pairing states are always unstable 
for s-wavc pairing. The resulting S-N phase diagram of 
this spin-singlet FS in the plane fi — h is presented in 
Fig. QJ, in which the dashed line indicates the boundary 
between the stable (left side) and unstable (right side) 
BCS states. 

The results for nonzero solutions of gap equation ([7]) 
for p+ip-wave symmetry are presented in Fig. 0^. At 
high fi the solutions have the same structure as in the 
s-wave case. They consist of Sarma states at high ex- 



change fields and BCS states with a critical exchange 
field h c — (fi/t)A p0 /V2. Decreasing fi causes a de- 
crease in the range of not only Sarma states but also 
BCS states. At a critical doping Sarma states are dis- 
appeared completely, and at the same time p-n solutions 
start to emerge. From this point on the gap equation 
support single- valued solutions consisting of BCS and p- 
n solutions, at low and high exchange fields, respectively. 
Upon approaching Dirac point there remains only single- 
valued p-n solutions, which extends up to a large critical 
field h c = A(l - 97Vi/2A)i, where Ni = 2Trv 2 F t 2 / gi and 
A is the highest energy scale as the cutoff in the integra- 
tion over the energy. As in the case of s-wave pairing, we 
have examined the stability of these solutions by plot- 
ting fls-N as a function of h for different values of /i. 
The results presented in Fig. 0b, show that in addition 
to BCS state whose range of stability (like the range of 
their existence) decreases with lowering fi, p-n states are 
also stable for all values of /i. We have presented the re- 
sulting phase diagram in Fig.QJl . As in Fig. [TJ;, we have 
depicted the boundary between the stable and unstable 
S states, whose dashed and solid parts indicate the first- 
and the second-order S-N phase transitions, respectively. 
We note that the above presented dependence of A SjP 
versus h can be used to distinguish between three differ- 
ent type of phases. Crossing a boundary between BCS 
and Sarma phases, the derivative of A SiP versus h will 
undergo a jump. While a p-n phase is distinguishable 
by its existence at exchange fields larger than the pairing 
gap, the fields for which a normal phase is expected in 
the ordinary FS. 

In conclusion, we have developed a mean field the- 
ory of relativistic-likc ferromagnetic superconductivity in 
graphene with an on-site and a nearest neighbor off-site 
pairing attraction. The resulted zero temperature phase 
diagrams for s-wave and p+ip-wave FS explore the ex- 
istence of a new gaplcss homogeneous S phase, which 
is formed by pairing of chiral carriers of opposite types, 
electron-like and hole-like. We have found that the p-n 
pairing with p+ip-wave symmetry presents a stable con- 
densate which, at low levels of doping, dominates the 
phase diagram up to exchange fields much larger than 
the superconducting pairing gap. 
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